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C*~> . The statistical objects characterizing turbulence in real turbulent flows differ 
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Q\ , tions from various two-dimensional and three-dimensional aspects, and from 



from those of the ideal homogeneous isotropic model. They contain contribu- 



the superposition of inhomogeneous and anisotropic contributions. We em- 

1. , ploy the recently introduced decomposition of statistical tensor objects into 

<S ■ 
r£^ irreducible representations of the S0(3) symmetry group (characterized by 

O 

j and m indices, where j = . . . oo, —j<m< j) to disentangle some of 

these contributions, separating the universal and the asymptotic from the 
specific aspects of the flow. The different j contributions transform differ- 
ently under rotations and so form a complete basis in which to represent the 
tensor objects under study. The experimental data are recorded with hot- 
wire probes placed at various heights in the atmospheric surface layer. Time 
series data from single probes and from pairs of probes are analyzed to com- 
pute the amplitudes and exponents of different contributions to the second 
order statistical objects characterized by j = 0, j = 1 and j = 2. The anal- 
ysis shows the need to make a careful distinction between long-lived quasi 
two-dimensional turbulent motions (close to the ground) and relatively short- 

1 



lived three-dimensional motions. We demonstrate that the leading scaling 
exponents in the three leading sectors (j = 0, 1, 2) appear to be different but 
universal, independent of the positions of the probe, the tensorial component 
considered, and the large scale properties. The measured values of the scaling 
exponent are C? =0) = 0.68±0.01, C? =1) = 1-0±0.15 and $ =2) = 1.38±0.10. 
We present theoretical arguments for the values of these exponents using the 
Clebsch representation of the Euler equations; neglecting anomalous correc- 
tions, the values obtained are 2/3, 1 and 4/3 respectively. Some enigmas and 
questions for the future are sketched. 
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I. INTRODUCTION 

The atmospheric boundary layer is a natural laboratory of turbulence that is unique in 
that it offers very high Reynolds numbers (Re). Especially if the measurements are made 
during periods when mean wind speed and direction are roughly constant, one approaches 
"controlled" conditions that are the goals of an experiment. Students of turbulence interested 
in the scaling properties, expected to be universal in the limit Re— > oo, are thus attracted 
to atmospheric measurements. On the other hand the boundary layer suffers inherently 
from strong inhomogeneity (explicit dependence of the turbulence statistics on the height) 
which leads to strong anisotropies such that the vertical and horizontal directions are quite 
distinguishable. In addition, one may expect the boundary layer near the ground to exhibit 
large-scale quasi two-dimensional eddies whose typical decay times and statistics may differ 
significantly from the generic three-dimensional motion. The aim of this paper is to offer 
systematic methods of analysis to resolve such difficulties, leading to a useful extraction of 
the universal, three-dimensional aspects of turbulence. 

Fundamentally, we propose to anchor the analysis of the statistical objects that are 
important in turbulence to the irreducible representations of the SO (3) symmetry group. 
Although the turbulence that we study is non-isotropic, the Navier-Stokes equations are 
invariant to all rotations. Together with incompressibility, this invariance implies that the 
hierarchy of dynamical equations satisfied by the correlation or structure functions are also 
isotropic |lj. This symmetry has been used in Q to show that every component of the 
general solution with a given index j, and hence a definite behavior under rotation, has 
to satisfy these equations individually, independent of other components with different be- 
havior under rotation. This "foliation" of the hierarchical equations motivates us to expect 
different scaling exponents for each component belonging to a particular j sector of the 
SO (3) decomposition. A preliminary test of these ideas has been reported in 0. The main 
result of the analysis shown below, supporting the results in , is that in each sector of the 
symmetry group scaling behavior can be found with apparently universal scaling exponents. 



We demonstrate below that scale-dependent correlation functions and structure functions 
can be usefully represented as a sum of contributions with increasing index j characterizing 
the irreducible representations of SO (3). In such a sum the coefficients are not universal, but 
the scale dependence is characterized by universal exponents. That is, a general component 
S a/3 (R) of the second rank structure function, 

S a0 (R) = ([u a (x + R)- u a (x)}[u^(x + R)- vP{x)\) , (1) 

where (. . .) stands for an ensemble average, can be usefully presented as a sum 

S<#(R) = £ a qjm {R)B a q U R ) , (2) 

where B q j m are the basis functions of the SO (3) symmetry group that depend on the direc- 
tion of R and a q j m are the coefficients that depend on the magnitude of R and, in general, 
include any scaling behavior. In other words, in terms of a rotation operator 0\ which 
rotates space by an arbitrary angle A one writes 

OkB$JR) = £ D^ m {K)B a J m {R) . (3) 

m'=-j 

The (2j + 1) x (2j + 1) matrices D~J m are the irreducible representations of the SO (3) 
symmetry group. The index q is required because, in general, there may be more than one 
independent basis function with the same indices j and m. We note that the basis functions 
Bqjm depend on the unit vector R only, whereas the amplitude coefficients a q j m (R) depend 
on the magnitude of R only. Our main point is that amplitudes scale in the inertial range, 
exhibiting universal exponents, 

a qjm (R) oc R^ . (4) 

Analyzed by usual log-log plots, a superposition such as (0) may well result in continuously 
changing slopes, as if there is no scaling. One of our main aims is to stress that the scaling 
exists, but needs to be revealed by unfolding the various contributions. This approach 
flushes out the scaling behavior even when the Reynolds number Re is low, as in numerical 
simulations ||. 



Obviously, to isolate tensorial components belonging to sectors other than the isotropic 
one needs to collect data from more than one probe. In Section 2 we present the experimental 
configuration and the conditions of measurement, and discuss the nature of the data sets. 
We demonstrate there that having two probes is actually sufficient to read surprisingly rich 
information about anisotropic turbulence. We have so far used two types of geometries, 
one consisting of two probes at the same height above the ground and the other with two 
probes vertically separated. In both cases the inter-probe separation is orthogonal to the 
mean wind. The caveat is that we must rely on Taylor's hypothesis Q to generate scale- 
dependent structure functions. In anisotropic flows the validity and optimal use of this 
method require discussion, and this is done in Section 3. In that Section we also examine 
the issues concerning 2- and 3-dimensional aspects of the flow pattern, and determine the 
outer scale L at which three-dimensional scaling behavior ceases to exist. In Section 4 we 
present the main results of the analysis. We demonstrate that the second-order structure 
function is best described as a superposition of contributions belonging to different sectors 
of the SO (3) symmetry group by extracting the coefficients and exponents that appear in 
the superposition (0). The scaling exponents depend on j, and we will demonstrate that 
they are an increasing function of j. For j = 0,1,2 our data analysis leads to the exponent 
values 0.68±0.01, 1.0±0.15 and 1.38±0.10 respectively. In Section 5 we present theoretical 
considerations that determine these exponents neglecting intermittency corrections, based 
on the Clebsch representation of the Euler equation. Section 6 offers a summary of the 
principal conclusions, and a discussion of the road ahead. 

II. THE EXPERIMENTAL SETUP 

The results presented in this paper are based on two experimental setups, which are 
denoted throughout as I and II respectively. In both setups the data were acquired over the 
salt flats in Utah with a long fetch. The site of measurements was chosen to provide steady 
wind conditions. The surface of the desert was very smooth and even. The measurements 
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FIG. 1. Schematic of the experimental set-up. Shown is the positioning of the probes with 
respect to the mean wind, and an explanation of how Taylor's hypothesis is employed. 

were made between 6 PM and 9 PM in a summer season during which nearly neutral stability 
conditions prevailed. The boundary layer was very similar to that on a smooth flat plate. 
In set I the data were acquired simultaneously from two single hot-wire probes at a height 
of 6 m above the ground, with a horizontal separation of 55 cm, nominally orthogonal to 
the mean wind, see Fig. [I]. 

The Taylor microscale Reynolds number was about 10,000. Set II was acquired from 
an array of three cross-wires, arranged above each other at heights 11 cm, 27 cm and 54 
cm respectively. The Taylor microscale Reynolds numbers for this set were 900, 1400 and 
2100 respectively. The hot-wires, about 0.7 mm in length and 6 /zm in diameter, were 
calibrated just prior to mounting them on the mounting posts and checked immediately 
after dismounting. The hot-wires were operated on DISA 55M01 constant-temperature 
anemometers. The frequency response of the hot-wires was typically good up to 20 kHz. The 
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TABLE I. Data sets I (first line) and II (second-fourth lines). 



voltages from the anemometers were suitably low-pass filtered and digitized. The voltages 
were constantly monitored on an oscilloscope to ensure that they did not exceed the digitizer 
limits. The Kolmogorov scales were about 0.75 mm (set I) and 0.5-0.7 mm (set II). Table I 
lists a few relevant facts about the data records analyzed here. The various symbols have the 
following meanings: U = local mean velocity, u' = root-mean-square velocity, (e) = energy 
dissipation rate obtained by the assumption of local isotropy and Taylor's hypothesis, r] 
and A are the Kolmogorov scale and Taylor microscale respectively, the microscale Reynolds 
number R\ = u'X/u, and f s is the sampling frequency. For set I we need to test whether 
the separation between the two probes is indeed orthogonal to the mean wind. (We do 
not need to worry about this point in set II, since the probes are above each other). To 
do so we computed the cross-correlation function (u\(t + r)u2(t)). Here, U\ and U2 refer to 
velocity fluctuations in the direction of the mean wind, for probes 1 and 2 respectively. If the 
separation were precisely orthogonal to the mean wind, this quantity should be maximum 
for t = 0. Instead, for set I, we found the maximum shifted slightly to r = 0.022 s, implying 
that the separation was not precisely orthogonal to the mean wind. To correct for this effect, 
the data from the second probe were time-shifted by 0.022 s. This amounts to a change 
in the actual value of the orthogonal distance. We computed this effective distance to be 
A t=a 54 cm (instead of the 55 cm that was set physically). We choose coordinates such that 
the mean wind direction is along the 3-axis, the vertical is in along the 1-axis and the third 
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direction orthogonal to these is the 2-axis. We denote these directions by the three unit 
vectors n, rh, and p respectively. The raw data available from set I is u^ 3 '(t) measured at 
the positions of the two probes. In set II each probe reads a linear combination of u^(t) 
and u^(t) from which each component is extractable. From these data we would like to 
compute the scale-dependent structure functions, using Taylor's hypothesis to surrogate 
space for time. This needs a careful discussion, which is given below. 

III. THEORETICAL CONSTRUCTS: TAYLOR'S HYPOTHESIS, INNER AND 

OUTER SCALES 

A. Taylor's Hypothesis 

Decades of research on the statistical aspects of hydrodynamic turbulence are based on 
Taylor's hypothesis [4-7], which asserts that the fluctuating velocity field measured by a 
given probe as a function of time, u(t), is the same as the velocity u(R/U) where U is the 
mean velocity and R = — Ut is the distance to a position "upstream" where the velocity is 
measured at t = 0. The natural limitation on Taylor's hypothesis is provided by the typical 
decay time of fluctuations of scale R. Within the classical scaling theory of Kolmogorov, 



this time scale is the turn-over time R/JS(R) where S(R) = S aa (R). With this estimate, 



Taylor's hypothesis is expected to be valid when JS(R)/U — ► 0. Since S(R) — * when 
R — > 0, the hypothesis becomes exact in this limit. We will use this aspect to match the 
units while reading a distance from a combination of space and time intervals. 

Reference presents a detailed analysis of the consequences of Taylor's hypothesis on 
the basis of an exactly soluble model. It also proposes ways for minimizing the systematic 
errors introduced by the use of Taylor's hypothesis. In light of that analysis we will use an 
"effective" wind, U e s, for surrogating the time data. This velocity is a combination of the 
mean wind U and the root-mean-square u', 



U cS = vV + (6m') 2 > (5) 
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where b is a dimensionless parameter. Evidently, when we employ the Taylor hypothesis in 
log-log plots of structure functions using time series measured in a single probe, the value 
of the parameter b is irrelevant, because it merely changes the (arbitrary) units of length 
(i.e, yields an arbitrary intercept). When we mix real distance between two probes and 
surrogated distance according to Taylor's hypothesis, the parameter b becomes a unit fixer. 
The numerical value of this parameter is found in Q by the requirement that the surrogated 
and directly measured structure functions coincide in the limit R — > 0. When we do not have 
the necessary data we will use values of b suggested by the exactly soluble model treated in 
0. This value of b ~ 3. We have checked that the scaling exponents change by no more 
than a few percent upon changing b by 30%. Further, this choice can be justified a posteriori 
by the quality of the fit of to the predicted scaling functions. 

When we have two probes placed at different heights, the mean velocities and u' as 
measured by the probes do not coincide. In applying Taylor's hypothesis one needs to 
decide the most appropriate value of U e R. This question has been addressed in detail in 
Ref. |7| with the final conclusion that the choice depends on the velocity profile between the 
probe. In the case of linear shear the answer is 
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where the subscripts 1,2 refer to the two probes respectively. 

In all subsequent expressions, we will therefore denote separations by R, and invari- 
ably this will mean Taylor-surrogated time differences or a combination of real and Taylor- 
surrogated distances. The effective velocity will be (|5|) or (^J) depending on whether the 
probes are at the same height or at different heights. 

B. Inner and Outer Scales in the Atmospheric Boundary Layer 

In seeking scaling behavior one needs to find the inner and outer scales. Below the 
inner scale second order structure functions have an analytic dependence on the separation, 
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FIG. 2. Log-log plots of the longitudinal component of the second order structure function. 

Panel (a) is for data set I and panel (b) for data set II. 

S(R) ~ R 2 , and above the outer scale they should tend to a constant value. We show in 
Fig. the longitudinal structure functions 

S 33 (R) = ((u®{x + R)- u {3 \x)) 2 ) (7) 

computed from a single probe in set I and from the probe at 0.54 m in set II. We also 
consider in Fig. R] the transverse structure function 



,M 



S lL (R) = ((u (l >(x + R)- u (L >(x)Y) 



(8) 



computed from the probe at 0.54 m in set II, see Fig. |3|. The spatial scales are computed 
using the local mean wind in both cases since we do not expect the scaling exponent for the 
single-probe structure function to be affected by the choice of advection velocity. However, 
this choice does determine the value of R corresponding to a particular time scale, but 
we expect that any correction to the numerical value of R is small for a different choice 
of advection velocity, and not crucial for the qualitative statements that follow. In Fig. |2| 
we clearly see the R 2 behavior characterizing the transition from the dissipative to the 
inertial range. As is well-known ||, this behavior persists for about a half-decade above the 
"nominal" Kolmogorov length scale rj. There is a region of cross-over and then the isotropic 
scaling ~ i? 0,68 expected for small scales in the inertial range begins. We thus have no 
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FIG. 3. Log-log plot of the transverse component of the second order structure function com- 
puted from data set II. 

difficulty in identifying the inner scale, it is simply revealed as a natural crossover length in 
these data. 

Next, since we cannot expect to fit with a single power-law for larger scales and must 
include scaling contributions due to anisotropy 0, we need to estimate the largest scales 
that should be included in our fitting procedure using the SO (3) decomposition machinery. 
We expect that the contributions due to anisotropy will account for scaling behavior up to 
the outer scale of the three-dimensional flow patterns. The task now is to identify that scale. 
One approach is simply to use the scale where the longitudinal structure function tends to 
a constant, corresponding to the scale across which the velocity signal has decorrelated. 
It becomes immediately apparent that this is not a reasonable estimate of the large scale. 
Figure || shows that the longitudinal structure function stays correlated up to scales that are 
at least an order of magnitude larger than the height at which the measurement is made. 
On the other hand, the transverse structure function computed from the probe at 0.54 m, 
Fig. |3], ceases to exhibit scaling behavior at a scale that is of the order of the vertical distance 
of the probe from the ground. 

It appears that we are observing extremely flat "quasi-two dimensional" eddies that are 
correlated over very long distances in the horizontal direction but have a comparatively small 
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vertical velocity component. Accordingly the vertical velocity component is dominated by 
bona fide three-dimensional turbulence. Since we know that the presence of the boundary 
must limit the size of the largest three-dimensional structures, the height of the probe should 
be something of an upper bound on the largest three-dimensional flow patterns that can be 
detected in experiments. Thus the size of the largest three-dimensional structures is more 
accurately determined by the decorrelation length of the transverse structure function. The 
theory of scaling behavior in three-dimensional turbulence can usefully be applied to only 
those flow patterns that are essentially three-dimensional. The extended flat eddies must 
be described in terms of a separate theory, including notions of two-dimensional turbulence 
which has very different scaling properties |§ . This is not the ambition of the present work. 
Rather, in the following analyses, we choose our outer-scale L in the horizontal direction 
to be of the order of the decorrelation length of the transverse structure function (where 
available) or of the height of the probe. We will see below that these two are the same to 
within a factor of 2; taking L to be as twice the height of the probe is consistent with all our 
data. We use this estimate in our study of both transverse, longitudinal and mixed objects. 

IV. EXTRACTING THE UNIVERSAL EXPONENTS OF HIGHER J SECTORS 

In this section we describe a procedure for extracting the scaling exponents that appear 
in the superposition fl2|). Preliminary results on the scaling exponent Q 7- \ obtained under 
the assumption of cylindrical symmetry, were announced in ||. The analysis here is more 
complete, and takes into account the full tensorial structure. We show that taking into 
account the full broken symmetry is feasible, and the final results are essentially the same. 
Both sets of results are also in agreement with analysis of numerical simulations ||. The 
results concerning Q 7- are new. 

In order to extract a particular j contribution and the associated scaling exponent, one 
would ideally like to possess the statistics of the velocity at all points in a three-dimensional 
grid. One could then extract the j contribution of particular interest by multiplying the 
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full structure function by the appropriate B q> j >m and integrating over a sphere of radius 
R. Orthogonality of the basis functions ensures that only the j contribution survives the 
integration. One could then perform this procedure for various R and extract the scaling 
behavior. This method was adopted successfully in || using data from direct numerical 
simulations. The experimental data are limited to a few points in space, so the integration 
over the sphere is not possible. We are faced with a true superposition of contributions from 
various j sectors with no simple way of disentangling them. However we can do the next 
best thing and use the postulate that the scaling exponents form a hierarchy of increasing 
values for increasing j. This can be interpreted to mean that anisotropic effects appear to 
increase with increasing scale. Since we look for the lowest order anisotropic contributions 
in our analyses, we perform a two-stage procedure to separate the various sectors. First we 
look at the small scale region of the inertial range to determine the extent of the fit with 
a single (isotropic) exponent. We then seek to extend this range by including appropriate 
anisotropic tensor contributions, and obtain the additional scaling exponents using least- 
squares fitting procedure. The following two sections discuss the procedure for determining 
the j '• = 2 and j = 1 scaling contributions to second-order statistics. 

A. The j=2 component 

In the second order structure function defined already, viz., 

S a ?(R) = ((u a (x +R)- u a {x)){u\x +R)- u (x))), (9) 

the j = 2 component of the SO (3) symmetry group corresponds to the lowest order 
anisotropic contribution that is symmetric in the indices, and has even parity in R (due 
to homogeneity). Although the assumption of axisymmetry used in |2j seemed to be justi- 
fied from the excellent qualities of fits obtained, we attempt to fit the same data (set I) with 
the full tensor form for the j '• = 2 contribution. The derivation of the full j '• = 2 contribution 
to the symmetric, even parity, structure function appears in Appendix A. 
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To begin with, we seek the range over which the isotropic scaling exponent holds for 
data set I. We measure all separation distances in units of A = 0.54 m which is the distance 
between the probes. The lower bound to the inertial range in this set is estimated to begin 
at R/A « 0.2 (see discussion in Section III B). We then find the range of scales over which 
the structure function 

S 33 (R, = 0) = ((4 3) (x + R)- uf\x)) 2 } } (10) 

with the subscript 1 denoting one of the two probes, can be fitted with a single exponent; 
this then would indicates the limit of isotropic scaling. We find that for points in the range 
0.2 < R/A < 4.5 a least squares fitting procedure yields a best-fit value (2 = 0.68 ± 0.01 
(Fig. H (a)). Fig. ^ (b) shows the fit to the structure function computed from a single probe 
in set I with just the j '• = contribution. Above this range, we are unable to obtain a 
good fit to the data with just the isotropic exponent and Fig. [| (b) shows the peel-off from 
isotropic behavior above R/A ~ 4.5. We point out that for ranges higher than this, one 
can indeed able to find a "best-fit" exponent for the curve, but the value of the exponent 
rapidly decreases and the quality of the fit is compromised. 

To find the j = 2 anisotropic exponent we need to use data taken from both probes. To 
clarify the procedure, we show in Fig. [I] the geometry of set I. What is computed is actually 

S 33 (R,6) = ([u^iUeet + U^-u^iUest)} 2 ) (11) 



Here 6 = arctan(A/[/ ef jt^), t^ = R/U e g, and R = J A 2 + (Ucst^) 2 . U e g is defined by Eq. ([|) 
with 6 = 3. For simplicity we shall refer from now on to such quantities as 

S 33 (R,9) = ((u i3 \x + R)-u { 2 3 \x)) 2 ) . (12) 

Next, we fix the scaling exponent of the isotropic sector as 0.68 and find the j = 2 
anisotropic exponent that results from fitting to the full j '• = 2 tensor contribution. We fit 
the objects in Eqs. fllPP and (|TJ) to the sum of the j = (with scaling exponent (2 = 0.68) 
and the j '• = 2 contributions (see Appendix A) 
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FIG. 4. The structure function computed from the single probe data, set I. (a) shows the 
X 2 minimization by the best-fit value of the exponent in the isotropic sector £2 ~ 0.68 for the 
single-probe structure function in the range 0.2 < R/A < 4.5. (b) shows the fit using the best 
value of C2 obtained in (a), indicating the peel-off from isotropic behavior at the end of the fitted 
range. 
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We fit the experimentally generated functions to the above form using values of Q 
ranging from 0.5 to 3. Each iteration of the fitting procedure involves solving for the six 
unknown, non-universal coefficients. The best value of Q is the one that minimizes the 
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FIG. 5. The x minimization by the best-fit value of the exponent in the j = 2 anisotropic 
sector from the fit to both the 6 = and the ^-dependent structure functions in the ranges 
0.2 < R/A < 25 and 1 < R/A < 25 respectively. 
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FIG. 6. The structure functions computed ^from data set I and fit with the j = and full 

(2) 

j = 2 tensor contributions using the best fit values of exponents £2 = 0.68 and Q = 1.38. Panel 
(a) shows the fit to the single-probe (0 = 0) structure function in the range 0.2 < R/A < 25 and 
panel (b) shows the fit to the ^-dependent structure function in the range 1 < R/A < 25. 
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TABLE II. The scaling exponents and the 6 coefficients in units of (m/sec) 2 as determined 
from the nonlinear fit of Eq. 7 to data set I. 

X 2 for these fits; from Fig. [5] we obtain this to be 1.38 ± 0.15. The fits with this choice of 
exponent are displayed in Fig. ^. The corresponding values of the six fitted coefficients is 
given in Table II. The range of scales that are fitted to this expression is 0.2 < R/A < 25 
for the 9 = (single-probe) structure function and 1 < R/A < 25 for the 9 ^ (two-probe) 
structure function. We are unable to fit with Eq. fll~3l ) to larger scales - - that is, larger 
than about 12 meters -- without losing the quality of the fit in the small scales. This is 
roughly twice the height of the probe from the ground. Based on the discussion in Sect. II 
C, we should be in the regime of the largest scales where the three-dimensional theory would 
hold. Therefore this limit to the fitting range is consistent with our expectations for the 
maximum scale of three-dimensional turbulence. We conclude that the structure functions 
exhibits scaling behavior over the whole scaling range, but this important fact is missed 
if one does not consider a superposition of the j = and j = 2 contributions. We thus 

(2) 

conclude that the estimate for the j ; = 2 scaling exponent Q ~ 1-38. This same estimate 
was obtained in [0] using only the axisymmetric terms. The value of the coefficients a and b 
are again close in magnitude but opposite in sign — just as in giving a small contribution 
to S 33 (R,9 = 0). The non-axisymmetric contributions vanish in the case of 9 = 0. The 
contribution of these terms to the finite 9 function is relatively small because the angular 
dependence appears as sin 9 and sin 2 9, both of which are small for small 9 (large R); and 
hence previously we were able to obtain a good fit to just the axisymmetric contribution. 
Lastly, we note that the total number of free parameters in this fit is 7 (6 coefficients and 
1 exponent). This brings up the possibility of having "over-fit" the data. The relative 
"flatness" of the \ 2 function near its minimum in Fig. |^ could be indicative of the large 
number of free parameters in the fit. However, the value of the exponent is perfectly in 
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agreement with the analysis of numerical simulations || in which one can properly integrate 
the structure function against the basis functions, eliminating all contributions except that 
of the j = 2 sector. Furthermore, fits to the data in the vicinity of Q — 1-38 show enough 
divergence from experiment that we are satisfied about the genuineness of the x 2 result. 

B. Extracting the j = l component 

The homogeneous structure function defined in Eq. (|9[) is known from properties of 
symmetry and parity to possess no contribution from the j = 1 sector (see Appendix B, 
Section 2), the j ' = 2 sector being its lowest order anisotropic contributor. In order to isolate 
the scaling behavior of the j = 1 contribution in atmospheric shear flows we must either 
explicitly construct a new tensor object which will allow for such a contribution, or extract it 
from the structure function itself computed in the presence of inhomogeneity. In the former 
case, we construct the tensor 

T a/3 (R) = (u a (x + R)- u a {x)){u\x + R) + vP{x))). (14) 

This object vanishes both when a = j3 and when R is in the direction of homogeneity. 
^From data set II we can calculate this function for non-homogeneous scale-separations (in 
the shear direction). In general, this will exhibit mixed parity and symmetry; we cannot use 
the incompressibility condition to reduce our parameter space. Therefore, to minimize the 
final number of fitting parameters, we examine only the antisymmetric contribution. We 
derive the tensor contributions in the j '• = 1 sector for the antisymmetric case in Appendix 
B, Section 1, and use this to fit for the unknown j = 1 exponent. We describe the results 
of this analysis below. For completeness, we have derived the tensor contributions in the 
j = 1 sector for the symmetric case as well in Appendix B2. This can be used to find 
j = 1 exponent for the inhomogenous structure function which is symmetric but has mixed 
parity. We do not present the results of that analysis here essentially because they are 
consistent with those from the antisymmetric case. Returning now to consideration of the 
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FIG. 7. The x 2 minimization by the best-fit value of the exponent Q °f the j = 1 anisotropic 
sector from the fit to ^-dependent T 31 (R, 9) function in the range 1 < R/A < 2.2. 
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antisymmetric part of the tensor object defined in Eq. fll4|), 



viz.. 



T a P{R) = - — ^ — - — ^ = (u a {x) u P{x + R)) - (u p (x)u a (x + R)), (15) 

which will only have contributions from the antisymmetric j = 1 basis tensors. An additional 
useful property of this object is that, it does not have any contribution from the isotropic 
j = sector spanned by 5 af3 and R a R@. This allows us to isolate the j = 1 contribution and 
determine its scaling exponent Q starting from the smallest scales available. Using data 
(set II) from the probes at 0.27 m (probe 1) and at 0.11 m (probe 2) we calculate 

f 31 (R) = {uf ) {x)u l i\x + R)) - (uf\x + R)ui\x)) , (16) 

where again superscripts denote the velocity component and subscripts denote the probe by 
which this component is measured. The goal is to fit this experimental object to the tensor 
form derived in Appendix Bl, Eq. (B7), 

f 31 (R,6,(f) = 0) = -03,1,0^ sinfl + a 2 ,i,iii: C ^ + a^^R^ cos 9. (17) 

Figure |7| gives the x 2 minimization of the fit as a function of Q . We obtain the best 
value to be 1 ± 0.15 for the final fit. This is shown in Fig. |8|. The fit in Fig. |^ peels off at 
around R/A = 2. The values of the coefficients corresponding to the exponent Q — 1 are 
given in Table III. The maximum range of scales over which the fit works is of the order 
of the height of the probes from the ground, consistent with the considerations presented 
earlier. This value of the scaling exponent of the j = 1 sector is entirely new. Again, we 
have satisfied ourselves that a different value of the exponent yields a substantially poorer 
fit to the data. 

In Section 5 we will present theoretical considerations to show that the value Q — 1 is 
predicted by a version of classical dimensional analysis. The present findings significantly 
strengthen our proposition |]J that the scaling exponents in the various sectors (at least up 
to j = 2) are indeed universal. 
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^2 


°3,1,0 


02,1,1 


03,1,-1 


1±0.15 


0.0116 ±0.001 


0.0124 ± 0.001 


-0.0062 ± 0.001 



TABLE III. The values of the exponents and coefficients (in units of (m/sec) 2 ) obtained from 
the fit to the function f 31 (R, 6). 

V. THEORETICAL DETERMINATION OF d J) FOR J = 1,2 



In this section we present dimensional considerations to determine the "classical K4T" 
values expected of Q and Q ■ We work at the same level as the K41 approach that yields 
the value Q = 2/3. This is justified since the differences between any two values Q an d 
Q J for j ^ j' are considerably larger than the intermittency corrections to either of them. 
We note, however, that the issue of anomalous exponents in turbulence has now multiplied 
several-fold, to all the j sectors, in light of the apparent universality that has unfolded in 
this work. 

It is easiest to produce a dimensional estimate for Q ■ One simply asserts || that the 
j = 2 contribution is the first one appearing in S al3 (R) due to the existence of a shear. 
Since the shear is a second rank tensor, it can appear linearly in the j '• = 2 contribution to 
S al3 (R). We thus write for any m, —j<m< j, 



Sf£ 2 (R)~T a W^f(R 7 e) 



(18) 



Here T a/37<5 is a constant dimensionless tensor made of 5 al3 , R a /R, and bilinear contributions 



made of the three unit vectors p, m, n, as exemplified in appendix A. The way Eq.(18) 
is represented means that the dimensional function f(R,e) stands for the response of the 
second order structure function to a small external shear. Ad such it is an inherent property 
of isotropic turbulence. Within the standard Kolomogorov-41 dimensional reasoning this 
function in the inertial interval can be made only of the mean energy flux per unit time and 
mass, e and R itself. The only combination of e and R that yields the right dimensions of 
the function / is e 1//3 i? 4 / 3 . Therefore 
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8TP 



S£f 2 (#) ~ T^^^e 1 / 3 /? 4 / 3 . (19) 

We thus find a "classical K41" value of Q = 4/3. Thus this simple argument seems to 

(2) 

rationalize nicely the experimentally found value Q = 1-38 ± 0.1. 

To understand the value of Q we cannot proceed in the same way. We need a contribu- 
tion that is linear (rather than bilinear) in the unit vectors p, m, n. We cannot construct a 
contribution that is linear in the shear and yet does not vanish due to the incompressibility 
constraint. Thus there is a fundamental difference between the j '• = 2 contribution and the 
j = 1 term. While the former can be understood as an inhomogeneous term linear in the 
forced shear, the j '• = 1 term, being more subtle, may be connected to a solution of some ho- 
mogeneous equation well within the inertial interval. In fact, all the known inertial-interval 
spectra in turbulent systems are related to the existence of a a flux of some conserved quan- 
tity which has a representation as an integral of some density in fc-space. For example the 
kinetic energy may be written as J dk\u(k, t)\ 2 . A well-known other integral of motion in 
hydrodynamics with such a presentations is the helicity 

H = f dr(u-V xu) . (20) 

Thus the helicity may be considered as a natural candidate which is responsible for a new 
solution in the inertial interval that may rationalize the j '• = 1 finding. We show that this is 
not the case in the following way. 

The dimensionality of H (denoted as [H]) differs from the dimensionality of the energy 
E by one length: [H] = [E/R]. Correspondingly, the dimensionality of the helicity flux, h 
may be written as 

[h] = [e/R] . (21) 

It means that in turbulence with energy and helicity fluxes one has at one's disposal a 
dimensionless factor in the form hR/e. It means that the second order structure function 
S(R, e, h) cannot be found just by dimensional reasoning even within the K41 approach. 
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Nevertheless, assuming that at small helicity fluxes (i.e. when hR/e <C 1) the function 
S(R, e, h) may be expanded in powers of h we can justify the first order correction due to 
helicity, 5hS, as 

S h S(R, e, h) ~ [e)- 1/3 hR 5/3 . (22) 

The value of the inferred scaling exponent, i.e. 5/3, is much larger than the value unity found 
experimentally. We thus need to find another invariant that may rationalize the findings. 

The only invariance in addition to the conservation of helicity that we are aware of in 
the inviscid limit is the Kelvin circulation theorem, which however does not furnish a local 
integral in fc-space in the Eulerian representation. The only way that is apparent to us to 
expose this invariance in a useful way is the Clebsch representation, in which one writes the 
Euler equation in terms of one complex field a(r, t), see for example |10|| . In fc-representation 



the Fourier component of the velocity field u(k, t) is determined from a bilinear combination 
of the complex field: 

u(k, t) = — [ d 3 M 3 &2*(fci, k 2 )a*(k u t)a(k 2 , t) , (23) 

07T d J 

*(k t , k 2 ) = Uk 1 + k 2 - fa - k 2 ) *'~^ ) . (24) 

It is well-known [[TIJ that this representation exposes a local conserved integral of motion 
which is 

n = — f d 3 kka*(k,t)a(k,t) . (25) 

on 6 J 

Note that this conserved quantity is a vector, and it cannot have a finite mean in an isotropic 
system. Consider now correction S n S(R, e, it) to the second order structure function due 
to a flux 7T of the integral of motion it. The dimensionality of it, [it] = [e 2 / 3 /R 1 ^ 3 ] and 
therefore now the dimesionless factor is it R 1 / 3 /e 2 / 3 . Assuming again expandability of S„S 
at small values of the flux it one finds that 

S^i{R) ~ T a ^ tt 7 R , (26) 
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where T a ^ is a constant dimensionless tensor linear in the unit vectors p, rh, h. We thus 
find the "classical K41" value Q = 1, which should be compared with the experimental 
finding C2 1} = 1 ±0.15. 

We stress that Eqs. ( pop and (|i~9| ) are the analogs of the standard isotropic dimensional 
estimate 

S?« (R) ~ (eR) 2/3 . (27) 

We thus conclude that dimensional analysis predicts that values 2/3, 1 and 4/3 for Q 
with j = 0, 1 and 2 respectively. This appears to be in satisfactory agreement with the 
experimentally extracted values of these exponents. We should state however that we do 
not know at present how to continue this line of argument for j > 2. 

VI. SUMMARY, CONCLUSIONS, AND THE ROAD AHEAD 

In summary, we considered the second order tensor functions of velocity in the atmo- 
spheric boundary layers. The following conclusions appear important: 

1. The atmospheric boundary layer exhibits three-dimensional statistical turbulence in- 
termingled with activities whose statistics are quite different. The latter are eddies 
with quasi two-dimensional nature, correlated over extremely large distances com- 
pared to the height of the measurement, having little to do with the three-dimensional 
fluctuations discussed above. 

2. We found that the "outer scale of turbulence" as measured by the three-dimensional 
statistics is of the order of twice the height of the probe. 

3. The inner scale is the the usual dissipative crossover, which is clearly seen as the scale 
connecting two different slopes in log-log plots. 

4. Between the inner and the outer scales, Eq. (|2|) appears to offer an excellent represen- 
tation of the structure function. Using contributions with j = 0, 1, 2 we could fit the 
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whole range very accurately. 

5. The scaling exponents Q J are measured as 0.68±0. 01, 1±0.15, 1.38±0. 10 respectively. 

6. Classical K41 dimensional considerations yield the numbers 2/3, 1 and 4/3 respectively. 
To get Q = 4/3 all that we need is to assume a contribution linear in the shear. For 
getting Q = 1 we need to identify a non-obvious conserved quantity which allows a 
new solution in the depth of the inertial interval. This is the first time that Clebsch 
variables allowed the understanding of a fundamentally new universal scaling exponent. 

If the trends seen here continue for higher j values, we can rationalize the apparent 
tendency towards isotropy with decreasing scales. If indeed every anisotropic contribution 

(i) 

introduced by the large scale forcing (or boundary conditions) decays as (R/L)^ with in- 
creasing Q as a function of j, then obviously when R/L — ► only the isotropic contribution 
survives. This is a pleasing notion that justifies the modeling of turbulence as isotropic at 
the small scales. 

We need to raise a word of caution here. We really have no idea about the values of 
the exponents for j > 3. Moreover, we are not even sure that they are well defined. To 
understand the difficulty one needs to examine the hierarchical equations for the correlation 
functions. These equations contain integrals used to eliminate the pressure contributions. 
The integrals were proven to converge (in the IR and the UV limits) when the the exponents 
(2 lies within the "window of convergence" which is (0, 4/3) (reference?). We see that with 
j = 3 we may reach beyond this window of convergence (this being questionable even for 
our experimental finding of j = 2!), and we are not guaranteed to have the kind of local 
theory that is thought to be a prerequisite to scaling behavior. 

Another enigma is related to the apparent success of the considerations of Section 5 
to rationalize the numerical values of the exponents found in the experiment. There is, 
however, no well-defined procedure of continuing the estimates for Q f° r J ' > 3. Whether 
this is related to the locality issue is not understood at present. 
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In conclusion, it appears that we have here an exciting possibility of generalizing the 
scaling structure of the statistical turbulence to many sectors of the symmetry group, gaining 
much better understanding of the structure of a theory. There exist, however, large patches 
of terra incognita on our map, patches that we hope to penetrate in future work. 
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APPENDIX A: FULL FORM FOR THE J = 2 CONTRIBUTION FOR THE 

HOMOGENEOUS CASE 

Each index j in the SO (3) decomposition of an n-rank tensor labels a 2 j + 1 dimensional 
SO(3) representation. Each dimension is labeled by m= — j, —j + 1, . . . j. The j = sector 
is the isotropic contribution while higher order j's should describe any anisotropy. The j = 
terms are well-known 



S*f (R) = c R^ 



(2 + ( 2 )5 ap - C 2 



(Al) 



R 2 

where (2 ~ 0.68 is the known universal scaling exponent for the isotropic contribution and cq 
is an unknown coefficient that depends on the boundary conditions of the flow. For the j '• = 2 
sector which is the lowest contribution to anisotropy to the homogeneous structure function, 
the m = (axisymmetric) terms were derived from constraints of symmetry, even parity 
(because of homogeneity) and incompressibility on the second order structure function M 
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S %, m ^R) = aR^ [(d 2) - 2)^ - C?\C? ] + 6) 
,(n^ (2) (2) OT^n-i?) 2 

+ ([d 2) ] 2 + 3d 2) + 6)nV 

a(2)/a(2)_ 2 n 

- ^ ^-g } -{R a n p + R^n a )(n-R)\ (A2) 

+ W* W [ - (d 2) + 3)(d 2) + 2)^(n -R? + ^- 

+ (d 2) + 3)(d 2) + 2)nV + (2d 2) + l)(d 2) - 2) 

x ^"^ • fi ) 2 _ ([C (2)]2 _ 4) (jR a n /? + H V)(n . #) 

(2) 

where Q is the universal scaling exponent for the j = 2 anisotropic sector and a and 6 are 
independent unknown coefficients to be determined by the boundary conditions. We would 
now like to derive the remaining m = ±1, and m = ±2 components 

Sj=2, m = 22 a q,2,mR 2 B"^ m (K) (A3) 

9 

where Q 3 ~ * s the scaling exponent of the j = 2 SO (3) representation of the n = 2 rank 
correlation function. The Bgj m (R) are the basis functions in the SO(3) representation of 
the structure function, The q label denotes the different possible ways of arriving at the the 
same j and runs over all such terms with the same parity and symmetry (a consequence 
of homogeneity and hence the constraint of incompressibility) [|TJ. In our case, even parity 
and symmetric in the two indices. In all that follows, we work closely with the procedure 
outlined in [|IJ. Following the convention in JTJ the o's to sum over are q = {1, 7, 9, 5}. The 
incompressibility condition d a u a = coupled with homogeneity can be used to give relations 
between the a q ,j,m for a given (j, m). That is, for j = 2, m = —2 ... 2 

(d 2) - 2)ai )2jr „ + 2(d 2) - 2)a 7Xm + (d 2) + 2)a 9 , 2 , m = (A4) 

0-1,2,171 + (C2 + 3)a7,2, m + C2 a 5,2,m = 

We solve equations (3) in order to obtain a 5i2 , m and a 7 ^, m m terms of linear combinations 
of Oi, 2 ,m and a g , 2 , m . 
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_ "-l,2,mUS2 J _ S2 ~ ^; "T "-9,2,mV[S2 J T ^S2 J_ ^ J ( a r \ 

a 5,2,m — ~ J .(2)/>(2) ~ l^^J 



a 7,2,r 



ai,2,m([d 2) ] 2 - d 2 ) -2) + g 9 , 2 , m ( [d 2) ] 2 + 5d 2) + 6) 

2d 2) (d 2) -2) 

Ol,2,m(2 — C2 ) — 09,2,m(2 + C2 ) 



2(d 2) - 2) 

Using the above constraints on the coefficients, we are now left with a linear combination 
of just two linearly independent tensor forms for each m 

Sj=2,m = a 9,2,m-R C2 [-(2 ( 2 + C2 )- B 7,2,m(- R ) + 2 C2 (C2 ~ 2 ) J B^2,m(- R ) 

+ ([d 2) ] 2 + 5d 2) + 6)5 5 t m ( J R)] 

+ a li2 , m ^ 2) [2d 2) (d 2) - 2)Bg >m (A) - Ci 2) (d 2) - 2)Bg m (A) 

+ ([d 2) ] 2 - d 2) - 2)B 5 lm(H)] (A6) 

The task remains to find the explicit form of the basis tensor functions Bg2 m (R), q G 
{l,7,9,5},me{±l,±2} 

• B$ im (R) = R- 2 8^WY 2m {R) 

• B$ im (R) = R- 2 [R a d? + R^d a ]R 2 Y 2m (R) 

• B^ 2 ,m(R) — R" A R a R^R 2 Y 2m (R) 

• Bti m {R) = d^d^R 2 Y jm {R) 

We obtain the m = {±1, ±2} basis functions in the following derivation. We first note 
that it is more convenient to form a real basis from the R 2 Y2 m (R) since we ultimately wish 
to fit to real quantities and extract real best-fit parameters. We therefore form the R 2 Y2k{R) 
(k = — 1, 0, 1) as follows: 

R 2 Y 2 o(R) = R 2 Y 2 o(R) = R 2 cos 2 9 = R 3 2 
R 2 Y 2 - 1 {R) = R 
= R 



*v ,£»_ d2 y 2 - 1 (R)-y 2+1 (R) 



2 
2 (cos — i sin 0) cos 9 sin 9 + (cos + i sin 0) cos 6* sin 



2 
= R 2 cos 6? sin 9 cos = -R3-R1 



/? 



-2i 

^cos 0— > % sin 0) cos 9 sin 6? — (cos + i sin 0) cos # sin 

^2i 
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= R cos 9 sin 9 sin = R%R 2 



2i 

2 (cos 20 + i sin 20) sin 2 — (cos 20 — i sin 20) sin 2 9 

= R % 

= R 2 sin 20 sin 2 9 = 2R t R 2 

2 (cos 20 + 2 sin 20) sin 2 6* + (cos 20 — i sin 20) sin 2 9 

= R 2 

= i? 2 cos 20 sin 2 9 = R\-R\ ( A7) 

This new basis of R 2 Y 2 k{R) is equivalent to using the R 2 Yj m (R) themselves as they form 
a complete, orthogonal (in the new k's) set. We omit the normalization constants for the 
spherical harmonics for notational convenience. The subscripts on R denote its components 
along the 1 (m) , 2 (p) and 3 (n) directions, m denotes the shear direction, p the horizontal 
direction parallel to the boundary and orthogonal to the mean wind direction and n the 
direction of the mean wind. This notation makes it simple to take the derivatives when we 
form the different basis tensors and the only thing to remember is that 

d a R 1 = d a (R • m) = m a 
d a R 2 = d a {R-p) =p a 
d a R 3 = d a (R-n)=n a (A8) 

We use the above identities to proceed to derive the basis tensor functions 

Bii-i( R ) = R~ 2 ^{R -n)(R- m) 

Bf^R) = R- 2 [{R a m p + R^m a ){R ■ n) + {R a n? + R^n a ){R ■ m)\ 

B$i-i(R) = R~ 2 R a R f3 {R -n)(R- m) 

B ti-ii.R) = nam/3 + n Q m a 

Bf 2A (R) = R~ 2 5 af3 {R -n)(R- p) 

Bf 2 1 (R)= R~ 2 [(R V + R V) (R-n) + (i? V + R p n a ) (R ■ p)] 
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BZiiiR) 
B?L 2 (R) 



R- 2 R a RP(R-n){R-p) 



n a pP _|_ n^p a 



B 



a/3 

7,2,-2 



} a/3 



B 



a/3 



'R) 



R) 



2R- 2 5 af3 {R-m){R-p) 

2R- 2 [(R a p /3 + R p p a )(R ■ m) + [R a mP + R p m a )(R ■ p)} 
BZM R ) = 2R~ 2 R a R p {R -m){R- p) 
= 2(mV + mV) 
= R- 2 5 al3 [(R-m) 2 -(R-p) 2 } 

= 2R~ 2 [{R a m p + R p m a )(R ■ m) - {R a p /3 + R p p a )(R ■ p)} 
= R~ 2 R a RP[(R-m) 2 -(R-p) 2 } 
= 2(m a m 13 - p a p (i ) 



5,2,-2 

Bf 2 , 2 (R) 



Bj^R) 
Bg ; 2,2(R) 



B t,2,2{R) 



(A9) 



Note that for each dimension k the tensor is bilinear in some combination of two basis 
vectors from the set m, p and n.bf. can we say something here about why this is so... using 
shear du a overdr f3 etc...?? Substituting these tensors forms into Eq. [A6| we obtain the full 
tensor forms for the j = 2 non-axisymmetric terms, with two independent coefficients for 
each k. 

' (2) r .(2)/ , r {2). 



S^ kss _ 1 (R) = OBA-iitf - Cf(2 + Cr)R- 2 [(R a m p + R?m a )(R 



n 



+ {R a n p + R p n a )(R ■ m)] + 2Q >(Q ' - 2)R~*R a R p (R ■ n)(R 



m) 



.(2) 



+ ([G f + 5Q } + 6)(n°W + n p m a ) 



>(2) 



+ 01,2,-1^ 2Cf } (Q ' - 2)R~ 2 6 a(i (R ■ n)(R 



m) 



*->j = 2,fc = U- ft ', 



- d 2) (d 2) - 2)R- 2 [{R a m (3 + R p m a ){R ■ n) + {R a n p + R p n a )(R ■ m)) 

+ ([d 2) ] 2 - d 2) - 2)(n a m 13 + n?m a ) 

= a 9 , 2 ,i^ 2) [ - d 2) (2 + d 2) )R- 2 [(R a pP + R p p a )(R ■ n) 

+ (R a n? + R?n a )(R ■ p)] + 2d 2) (d 2) - 2)R~ 4 R a R^(R -n)(R- p) 

+ ([d 2) ] 2 + 5d 2) + 6)(nV + nV)' 

+ a li2il i^ 2) [2d 2) (d 2) - 2)iT 2 <^(# • n)(R ■ p) 
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- d 2) (d 2) - 2)R- 2 [(R a pP + RPp a )(R ■ n) + (R a n p + R^n a )(R ■ p)] 
+ ([d 2) ] 2 - d 2) - 2)(nV + nV)' 



S&^iR) = 09,2,-2^ " - 2d 2) (2 + d 2) )i?- 2 [(i?V + i#V)(Jl 



TO 



+ (/Trr/ + R p m a )(R ■ p)} + 2Q J (Q } - 2)R~ i R a R p (R ■ p)(R 



TO 



+ ([d 2) ] 2 + 5d 2) + 6)(toV + mV)' 
+ oi A _ 2 i2^ [2d 2) (d 2) - 2)RT 2 S°*(R -m)(R- p) 
- 2d 2) (d 2) - 2)iT 2 [(#V + R P p a ){R ■ m) + (iFrn^ + R?m a ){R ■ p)] 
+ 2([d 2) ] 2 - d 2) - 2)(toV + mV)' 



S 



a/3 
j=2,k=2 



R 



09,2,2 



& 



(2) 



^•P) S 



2d 2) (2 + (i 2) )R- 2 l(R a mP + R^m a )(R ■ to) 

- (i?V + R P p a ){R ■ p)] + 2d 2) (d 2) - 2)R- i R a R f5 [(R ■ mf 
+ 2([d 2) ] 2 + 5d 2) + 6)(m a m fi - p p p a ) 

+ a 1)2 , 2 i^ 2) [2d 2) (d 2) - 2)iT 2 <^[(# • to) 2 - (# • pf] 

- 2d 2) (d 2) - 2)i?- 2 [(i? a TO /3 + R?m a ){R • to) - (i?V + R^p a )(R • p)] 



2(td 2) ] 2 



a(2) 
S2 



2)(mV-py) 



(A10) 



(2) 
Now we want to use this form to fit for the scaling exponent Q i n t he structure function 

S 33 (R) from data set I where a = /3 = 3 and the azimuthal angle of R in the geometry is 

= 7T/2. 

Sg 2 , fc= _ 1 (i2,e,0 = 7r/2) = 



42), 



S™ 2tk=1 (R, 6,(j> = tt/2) = 09,2,1^" [-2d 2) (d 2) + 2) sin cos 

+ 2d 2) (d 2) -2)cos 3 ^sin^] 



S£ 2 , fc= _ 2 (i2,M = */2) = 



(2), 



S 3 i 2>fc=2 (i?, 0, = tt/2) = 09,2,2^ [-2Cr ; (Q ; - 2) cos 2 sin 2 9} 



>(2) 



+ ai A2 i#'[-2Cncr ; -2)siii 2 0] 



(All) 



We see that choosing a particular geometry eliminates certain tensor contributions. In the 
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<j) = vr/2, 


a = j3 = 2> 


4> = 0, a = j3 = 2> 


= O,a = /3 = l 


= O,a = 3,/? = l 


6 + 


6> = 


6 + 


6 = 


fly^O 


6 = 


fly^O 


9 = 
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2 
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2 





-1 
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1 
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2 


2 





2 





2 


2 


2 





Total 


5 


2 


5 


2 


5 


4 


6 


2 



TABLE IV. The number of free coefficients in the j = 2 sector for homogeneous turbulence 
and for different geometries 

case of set I we are left with 3 independent coefficients for m^O, the 2 coefficients from the 
m = contribution (Eq. |A2| ), and the single coefficient from the isotropic sector |Al|, giving 



a total of 6 fit parameters. The general forms in [A10| can be used along with the k = 
(axisymmetric) contribution [Al] to fit to any second order tensor object. For convenience, 
Table IV shows the number of independent coefficients that a few different experimental 
geometries we have will allow in the j = 2 sector. It must be kept in mind that these forms 
are to be used only when there is known to be homogeneity. If there is inhomogeneity, then 
we cannot apply the incompressibility condition to provide constraints in the various parity 
and symmetry sectors and we must in general mix different parity objects, using only the 
geometry of the experiment itself to eliminate any terms. 



APPENDIX B: THE J=l COMPONENT IN THE INHOMOGENEOUS CASE 



1. Antisymmetric contribution 



We consider the tensor 



T af3 (R) =< u a (x + R)- u a {x)){u p {x + R)+ u\x)) > 



(Bl) 
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. This object is trivially zero for a — /3. In our experimental setup, we measure at points 
separated in the shear direction and therefore have inhomogeneity which makes the object 
of mixed parity and symmetry. We cannot apply the incompressibility condition in same 
parity/symmetry sectors as before to provide constraints. We must in general use all 7 
irreducible tensor forms. This would mean fitting for 7x3 = 21 independent coefficients 
plus 1 exponent Q m the anisotropic sector, together with 2 coefficients in the isotropic 
sector. In order to pare down the number of parameter we are fitting for, we look at the 
antisymmetric part of T a ^(R) 

T a(3 {R) = ^ ^ = (u a (x)u (x + R)) - (u p {x)u a {x + R)) (B2) 

which will only have contributions from the antisymmetric j '• = 1 basis tensors. These are 

• Antisymmetric, odd parity 

Bti m = R-^ITdP - R^d a }RY hm (R) (B3) 

• Antisymmetric, even parity 

^4,l,m = R e ° ^R^Ry\,m{R) 

Bf^m = R- 2 e a ^RY 1 , m (R) (B4) 

As with the j = 2 case we form a real basis RYi^(R) from the (in general) complex 
RYi^ m (R) in order to obtain real coefficients in our fits. 

RY ljk=0 (R) = RY lfi (R) = RcosO = R 3 

Y hl (R)+Y hl (R) 



RY hk=1 (R) = R 



RY 1 , = UR) = R 



2i 
R sin 9 sin <p = R 2 

Y 1 ,(R)-Y 1 AR) 



2 
R sin 9 cos <fi = R\ 



(B5) 



And the final forms are 
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B™%{R) = FC l [R a nP - R 13 



n~ 



Bf lfi (R)=R- 2 e a ^R fl (R. 



n) 



B%lo(R) = R~ 2 e a ^ 



n„ 



B^ l:1 (R) = R~ l [R a ^ - RPp a ] 

B2l 1 (R)=R-\^R,(R.p) 

Bf hl (R) = RT\^ ViX 
B^_ X {R) = R~ l \R a m p - R?m a ] 
B%l_ x (R) = R- 2 e a ^R^R.m) 
B^ 1 {R)=R- 2 e a ^m, (B6) 

Note: For a given k the representations is symmetric about a particular axis in our chosen 
coordinate system (l=m (shear), 2=p (horizontal), 3=n (mean-wind)) 

We now have 9 independent terms and we cannot apply incompressibility in order to re- 
duce the number of independent coefficients in our fitting procedure. We use the geometrical 
constraints of our experiment to do this. 

• = (vertical separation), a = 3, (3 = 3 

BH h0 (R,9,<p = 0) = -sm9 

B% 1A (R,e,<i> = 0) = 1 

B 3A,-i( R > M = 0) = cos 6 (B7) 

There are no contributions from the reflection-symmetric terms in the j = isotropic 
sector since these are symmetric in the indices. The helicity term in j = also doesn't 
contribute because of the geometry. So, to lowest order 

f a(3 (R)=f^ 1 (R) 

= a 3> i )0 (i2)(- sin 9) + 02,1,1 (i?) + 03,1,-1 (.R) cos 9 (B8) 

We have 3 unknown independent coefficients and f unknown exponent to fit for in our 
data. 
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2. Symmetric contribution 

We consider the structure function 

S af3 {R) =< (u a (x + R)- u a (x))(u^(x + R)- vP(x)) > (B9) 

in the case where we have homogeneous flow. This object is symmetric in the indices by 
construction, and it is easily seen that homogeneity implies even parity in R 

S al3 (R) = S Pa (R) 
S afS (-R) = S a(S (R) (BIO) 

We reason that this object cannot exhibit a j — 1 contribution from the SO (3) representation 
in the following manner. Homogeneity allows us to use the incompressibility condition 

d a S a/3 = 

dpS af3 = (Bll) 

separately on the basis tensors of a given parity and symmetry in order to give relationships 
between their coefficients. For the even parity, symmetric case we have for general j > 2 just 
two basis tensors and they must occur in some linear combination with incompressibility 
providing a constraint between the two coefficients. However, for j = 1 we only have 
one such tensor in the even parity, symmetric group. Therefore, by incompressibility, its 
coefficient must vanish. Consequently, we cannot have a j = 1 contribution for the even 
parity (homogeneous), symmetric structure function. Now, we consider the case as available 
in experiment when R has some component in the inhomogeneous direction. Now, it is 
no longer true that S al3 (R) is of even parity and moreover it is also not possible to use 
incompressibility as above to exclude the existence of a j = 1 contribution. We must look 
at all j '• = 1 basis tensors that are symmetric, but not confined to even parity. These are 
• Odd parity, symmetric 

Bf lk (R) = R- l 5 a ?RY lk (R) 
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B iik(R) = FT l [R a d p + R?d a ]RY lk (R) 

B^i k (R) = R^R a R p RY lk {R) 

Bf 1:k (R) = Rd a d?RY lk (R) = (B12) 

• Even parity, symmetric 

Bf 1>k (R) = R-2[R a e^R^d v + R^e a ^ R^d u ]RY lk (R) 

Beik(R) = l^ v R»d v d a + e^R,d v dp\RY lk (R) = (B13) 

We use the real basis of R^Yi^R) which are formed from the i? -1 Yi m (.R). Both B^ lk (R) 
and .Bgi k (R) vanish because of the taking of the double derivative of an object of single 
power in R. We thus have 4 different contributions to symmetric j '• = 1 and each of these is 
of 3 dimensions (k = — 1, 0, 1) giving in general 12 terms in all. 

Bf lfi (R) = R- 1 5 a/3 (R-n) 

B%%{R) = BT^IP-nP + flV] 

Bl%(R) = R~ 3 R a RP(R ■ n) 

Bsio(R) = R~ 2 [{R a mP + RP m a )(R ■ p) - (i?V + R^p a ){R ■ m)\ 

Bf 1:1 (R) = R- 1 5 a/3 (R-p) 

B$ A (R) = R- l [R a p? + i?V] 

Bf hl (R) = R-'R a R^R-p) 

B$ tl (R) = R- 2 [(R a m (3 + RP m a )(R ■ n) - (R a n fi + R p n a ){R ■ m)\ 
Bf lt _ 1 (R) = R- 1 5 a ^(R-m) 
Bf^R) = R-^m 13 + R?m a ] 
B^iR) = R- 3 R a R?(R ■ m) 
Bsi-i(R) = R' 2 [{R a p P + R p p a ){R ■ n) - (R a n P + RPn a )(R ■ p)\ (B14) 

These are all the possible j '• = 1 contributions to the symmetric, mixed parity (inhomo- 
geneous) structure function. 
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k 


(p = 0, a = (3 = 3 


(f) = 0,a = P = l 


= O,a = 3,/? = l 


0/0 


= 


0/0 


= 


0/0 


= 





3 


3 


2 


1 


2 





1 


1 





1 











-1 


2 





3 





2 


1 


Total 


6 


3 


6 


1 


4 


1 



TABLE V. The number of free coefficients in the symmetric j = 1 sector for inhomogeneous 
turbulence and for different geometries. 

For our experimental setup II, we want to analyze the inhomogeneous structure function 
in the case a — (3 — 3, and azimuthal angle = (which corresponds to vertical separation) 
and we obtain the basis tensors 



b\IM 


= cos 9 


S 7 3 l,o(#) 


= 2 cos e 


S 9 3 l,o(#) 


= cos 3 9 


BUM 


= — 2 cos # sin # 


m,-M 


= sin# 


m,-M 


= cos 2 9 sin 9 



(B15) 

Table V gives the number of free coefficients in the symmetric j '• = 1 sector in the fit to 
the inhomogeneous structure function for various geometrical configurations. 
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